Astronomy & Astrophysics manuscript no. final version-arxiv 


©ESO 2012 


December 19, 2012 





A new nearby pulsar wind nebula overlapping the RX J0852.0-4622 

supernova remnant 

F. Acero^'^, Y. Gallant\ J. Ballet^ M. Renaud^andR. Terrier"^ 

^ LUPM, CNRS/Universite Montpellier 2 e-mail: facero@in2p3 . fr 

2 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 

3 AIM, CEA/CNRS/Universite Paris Diderot-Paris 7 
^ APC, CNRS/Universite Paris Diderot-Paris 7 

Received 2012 November 27; accepted 2012 December 17 

ABSTRACT 

Context. Energetic pulsars can be embedded in a nebula of relativistic leptons which is powered by the dissipation of the rotational 
energy of the pulsar. The object PSR J0855-4644 is an energetic and fast-spinning pulsar (E = 1.1 x 10^^ erg/s, P=65 ms) discovered 
near the South-East rim of the supernova remnant (SNR) RX J0852.0-4622 (aka Vela Jr) by the Parkes multibeam survey. The posi- 
tion of the pulsar is in spatial coincidence with an enhancement in X-rays and TeV y-rays, which could be due to its putative pulsar 
wind nebula (PWN). 

Aims. The purpose of this study is to search for diffuse non-thermal X-ray emission around PSR J0855-4644 to test for the presence 
of a PWN and to estimate the distance to the pulsar. 

Methods. An X-ray observation was carried out with the XMM-Newton satellite to constrain the properties of the pulsar and its neb- 
ula. The absorption column density derived in X-rays from the pulsar and from different regions of the rim of the SNR was compared 
with the absorption derived from the atomic (HI) and molecular (^^CO) gas distribution along the corresponding lines of sight to 
estimate the distance to the pulsar and to the SNR. 

Results. The observation has revealed the X-ray counterpart of the pulsar together with surrounding extended emission thus confirm- 
ing the existence of a PWN. The comparison of column densities provided an upper limit to the distance of the pulsar PSR J0855-4644 
and the SNR RX J0852.0-4622 (d< 900 pc). Although both objects are at compatible distances, we rule out that the pulsar and the 
SNR are associated. With this revised distance, PSR J0855 -4644 is the second most energetic pulsar, after the Vela pulsar, within a 
radius of 1 kpc and could therefore contribute to the local cosmic-ray e~/e^ spectrum. 

Key words. Pulsars: general; ISM: supernova remnants; ISM: individual objects: (PSR J0855-4644, RX J0852.0-4622); 



1. Introduction 

The Vela region (approximately centered on ^ = 265°, b = -1°) 
is the host of several remnants of stellar explosions including 
the eponymous supernova remnant (SNR). The Vela remnant is 
among the brightest and most extended X-ray sources in the sky 
and it is also the closest SNR currently detected (the di stance to 
the remnant 's pulsar has been measured to be 290 pc; Dodson| 
|et al.[[2Q03] ). With its large angular size (8° in diameter) the ob- 
ject is overlapping several remnants such as the Puppis A SNR 
and RX J0852.0-4622. In the X-ray spectral domain, this com- 
plex region is dominated by the Vela and Puppis remnants in soft 
X-rays and the RX J0852. 0-4622 SNR (aka V ela Jr) becomes 
apparent only above 1 keV (| Aschenbach| 1 1 998| . 



The Parkes multibeam radio survey discovered a pulsar, 
PSRJ0855-4644, l ying on the South-E ast rim of the RX 
J0852.0-4622 SNR dKramer et aLj [2003] ). This young and en- 
ergetic pulsar (characteristic age Tc=140 kyrs, E = 1.1 x 10^^ 
erg/s and P=65 ms) is in spatial coincidenc e with an enhance 
ment in the X- and y-ray maps (presented in [Aschenbach et aF 



1999[|Aharonian et al.| 12007) respectively), which could be the 



signature of its pulsar wind nebula. Based on the r adio dispersion 



measure of the pulsar (DM=239 cm ^ pc; Kramer et al. 20031) 



combined with the NE2001 model of Galactic electron distribu- 
tion from Cordes & Lazio ( 2002| ), the distance to the pulsar is 



estimated to be 4 kpc. However, the estimation of the distance 
using such methods in this specific region of the sky can be erro- 
neous (see [Mitra & Ram achandran , 2001 ) due to the structured 
and nearby emission of the Gum nebula, a 35° large nebula lo- 
cated at -500 pc. The electro n density inside the Gum nebula can 
vary from 0.1 to 100 cm"^ (Reynolds 1976 IWallerstein et al. 



1980; Dubn er etaLlp^92l[Reynoso & Dubner , 1997 ) and there- 



fore so can the DM of a pulsar lying in the background of the 
nebula. This complex region is modeled as a simple sphere (the 
Gum nebula) of homogenous density in the NE2001 model. 

In order to better constrain the distance to the pulsar, and 
search for a potential X-ray nebula, an XMM-Newton observa- 
tion was carried out whose results are presented in this arti- 
cle. In Sect. |2) we present a detailed morphological and spec- 
tral analysis of the pulsar and of its nebula. A discussion of the 
distance of the PSR and of the SNR RX J0852.0-4622 is pre- 
sented in Sect. [3] To conclude, a summary of the proposed sce- 
nario, the potential relation between PSR J0855 -4644 and the 
RX J0852.0-4622 SNR and the expected TeV emission from 
the PWN is presented in Sect. [4] 

2. X-ray observations 

The SNR RX J0852.0-4622 was first discovered in X-rays by 
the ROS AT satellite in a high energy image (E> 1 .3 keV) and has 
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Fig. 1. XMM-Newton background subtracted flux images combining MOS and PN cameras in units of ph/cm^/s/arcmin^. Left: low energy band 
image (0.5-0.8 keV) adaptively smoothed to a signal-to-noise ratio of 10 using the task asmooth provided in the XMM-Newton software SAS. 
The position of PSR J0855-4644 is shown with a red circle. Right: image in the high energy band (1.2-6 keV) smoothed with a Gaussian width 
of 7". The dashed red line shows approximately the shell of the SNR RX J0852. 0-4622. In solid line we show the spectral extraction regions 
used to study the rim of RX J0852.0-4622 and the PWN whereas the dashed line regions are used to estimate the spectrum of the astrophysical 
background. 



revealed the energy dependent morphology of this complex re- 
gion. Due to its large angular size (nearly 1° in radius), the SNR 
has only been partly observed by the new generation of X-ray 
telescopes and, so far, no mosaic, covering the whole remnant, 
at high angular resolution is available. Most observations have 
been focused on the intriguing bright and thin filaments emitting 
synchrotron emission in the North- West of the remnant (see e.g. 
IBamba et al.||2QQ5] ). 

The XMM-Newton observation presented in this study is cen- 
tered on the pulsar and also unveils the, previously unobserved. 
South-eastern shell of the RX J0852.0-4622 SNR. Several other 
XMM-Newton observations in the North-West and in the South 
of the SNR were used in this study to compare the spectral prop- 
erties of the SNR along the shell. 



fuse emission is also seen surrounding the pulsar. The spectral 
properties of the pulsar and of the diff'use emission are presented 
in the following sections. 

The thermal emission from the Vela SNR is extremely bright 
at low energy (E<1 keV) and therefore the statistics gathered 
with the PN camera is very high. We have observed a calibration 
issue with the spectral response of the PN camera, particularly 
visible in the Oxygen lines. This is a known issue for the XMM- 
Newton team and following their recommendation an off'set in 
energy in the PN spectral response function {rmf) was intro- 
duced via the gain function in Xspec. All spectra were fitted us- 
ing MOS 1 +2 and PN. For graphical purposes, only the MOS 1 +2 
spectra are shown along with the best-fit obtained with the three 
instruments. See Appendix [A| for more details. 



2.1. Data analysis 

The region of the pulsar was observed with XMM-Newton for 55 
ks using the Full Frame mode for the three EPIC instruments. 
After flare screening, the remaining exposure time is 40 ks and 
29 ks for the MOS and PN cameras respectively. For the im- 
age generation, the instrumental background was derived from 
a compilation of blank sky observations ( jCarter & R ead , 2007) 
renormalized in the 10-12 keV energy band. The background 
subtracted and vignetting corrected images in a low (0.5-0.8 
keV) and high energy band (1.2-6 keV) are presented in Fig. 
[T] Those images clearly show two very difl'erent facets of the 
Vela region. While the low energy image is dominated by the 
large scale structure of thermal emission from the Vela SNR, the 
high energy band reveals a faint rim of the RX J0852.0-4622 
SNR and the X-ray counterpart of the pulsar. Interestingly, dif- 



2.2. Morphology of the nebula 

To further investigate the morphology of this difl'use emission, 
a radial profile (shown in Fig. 2j middle), centered on the posi- 
tion of the pulsar, was extracted in the high energy band (1.2-6 
keV). The profile was then fitted with a combination of a point- 
like component (using the XMM-Newton PSF from the calibra- 
tion files), an extended component (a Gaussian profile) and a 
constant respectively representing the pulsar, the nebula and the 
background. The best-fit Gaussian width deconvolved from the 
PSF of the extended component is 45"+ 1.5" (unless otherwise 
noted, all error bars are given at a 90% confidence level) and the 
diff'use emission around the pulsar is detected up to -150". The 
large scale emission is fairly symmetric and no indications of a 
trail or of a bow- shock structure (reminiscent of a fast moving 
pulsar) are observed. 
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Interestingly, at smaller scale the inner part of the nebula (r 
< 60") seems to show some azimuthal variation. To test the sig- 
nificance of these variations, an azimuthal profile was extracted 
from an annulus region (Rl = 15"< r <R2=60", see Fig. |2] top 
panel for region definition). The resulting brightness distribu- 
tion presented in Fig. [2] (bottom) is not compatible with a flat 
profile and exhibits features, 180° apart, at 50° and 230° which 
are suggestive of jet-like structures as observed in many other 
young pulsars (see |Kargaltsev & Pavlovl ^008 ). However, the 
limited spatial resolution of XMM-Newton hampers a study of 
the immediate environment of the pulsar. To confirm the pres- 
ence of jets in the nebula and to search for a torus, high spatial 
resolution observation with the Chandra X-ray satellite would 
be needed. 



2.3. Spectral analysis of the nebula 

To extract the spectral properties of the nebula, the instrumen- 
tal and astrophysical backgrounds have to be modeled carefully 
in this complicated region of the sky. The instrumental and par- 
ticle induced background were derived from a compilation of 
observations with the filter wheel in closed positiorQand renor- 
malized in the 10-12 keV energy band over the whole field of 
view. The spectral data were fitted using Cash statistics and the 
goodness-of-fit was assessed trough Monte Carlo simulations 
using the goodness command in Xspec (vl2.5). This command 
returns the fraction of simulated spectra whose fit statistics value 
is less than the one obtained on the real data which is equivalent 
to l-p-valu^ If the model is a good representation of the data, 
a p-value around 0.5 is expected. If the p-value value is close to 
0, then the model is clearly a bad representation of the data. The 
Tuebingen-Boulder ISM absorption model (tbabs) is used for all 
absorbed co mponents together with the abundances of [Wilms 
|etaLl ( [2QQQl ). 

The astrophysical background in this region is particularly 
complex as illustrated in Fig.[T] The low energy part of the spec- 
trum is dominated by the thermal emission from the Vela SNR 
and from the Local Bubble while in the high energy band non- 
thermal emission from RX J0852. 0-4622 (the contribution from 
the SNR's rim at the PWN's position is not negligible; see Sect. 
|2.6| ) and the cosmic X-ray background (CXB) are present. To 
model this background a spectrum from a region (dashed circle 
in Fig.[T] right) located at the same distance from the shock front 
of the RX J0852.0-4622 SNR (to have similar spectral proper- 
ties) has been obtained. 

The spectrum from this background region is well fitted by 
a two component model (p-value=0.24). The first component is 
an absorbed VAPEC (Variable abundances astrophysical plasma 
emission code) model which represents the thermal contribution 
from Vela. The X-ray absorption column A^h of this component 
is not well constrained with XMM-Newton' s energy range (0.4- 
10 keV) and is therefore fixed to 5 x 10^^ cm"^, the typical value 
derived by 
satellite w 



sorbed power-law model representing the contributions from RX 
J0852.0-4622 and the CXB. Adding another component to rep- 
resent separately the CXB, with spectral parameters fixed from 



Aschenbach et al. (1995) using data from the ROSAT 



^hose energy range extends to lower energy (0.1-2.4 
keV). The best-fit parameters for this thermal component are the 
following: kT=0.14 keV, Oxygen abundance=0.4^ and a nor- 
malization of 7.1 X 10"^ cm"^. The second component is an ab- 



^ http://xmm2.esac.esa jnt/external/xmm_sw_cal/background/filter_closed/ 
^ The p-value is the probability of obtaining a fit statistic greater than 

the observed value. 
^ The abundance of Neon was also let free as a test but was found 

to be compatible with solar abundance and was therefore fixed to solar 

value in the model. 



cm ^, 
^ keV"^cm"^ s"^ at 



De Luca & Molen di ( 2004), did not improve the quality of the 
fit. The best-fit parameters give A^h=(0.45 ± 0.10) x 10^^ 
r=1.90±0.04 and norm=(2.70 ± 0.11) x 10 ' 
1 keV. 

The best-fit model obtained from the background spectrum 
is then used as a fixed template model in the spectrum of the 
nebula (annulus: Rl = 15"< r <R3=150") where the normaliza- 
tion is set to the ratio of geometrical areas for the non-thermal 
emission and let free for the thermal emission as the emission 
from the Vela SNR is seen to slightly vary in the low energy 
map shown in Fig. [T] All the other parameters (kT, power-law 
absorption column and spectral index) for the background model 
are kept fixed. To represent the emission from the nebula, an ab- 
sorbed power-law model was added to the previously derived 
background model. The spectrum from the nebula is well mod- 
eled (see Fig.|3] top) with the sum of the background -h power- 
law model (p-value=0.14) thus confirming the non-thermal na- 
ture of the nebula's X-ray emission. 



2.4. Spectral analysis of PSRJ0855-4644 

In order to derive the X-ray properties of the pulsar, and in par- 
ticular the column density along the line of sight to discuss its 
distance, we have extracted a spectrum in a 15" radius region 
centered on the pulsar. For this region, the astrophysical back- 
ground is composed of Vela and the Local Bubble at low ener- 
gies and of the sum of the contributions from the nebula (which 
is not negligible based on the radial profile modeling in Fig. [2]), 
RX J0852.0-4622 and the CXB. 

The astrophysical background model was obtained from the 
best-fit model of the spectrum extracted in the PWN region 
(15"< r <150"). The non-thermal component of the template 
was renormalized using the nebula's geometrical model shown 
in Fig. [2] (middle) while the thermal component was renormal- 
ized using the geometrical ratio of the areas (i.e. thermal emis- 
sion is assumed to be homogeneous at the scale of the PWN). 

The spectrum of the pulsar (presented in Fig. |3j Bottom) is 
well described (p-value=0.46) at low energies by the fixed tem- 
plate (absorbed VAPEC model) and at high energies with two 
absorbed power-laws one representing the astrophysical back- 
ground and one for the non-thermal emission of the pulsar. The 
spectral parameters of the pulsar best-fit model are listed in Table 
[T] A black body model has also been tested to model the X-ray 
emission of the pulsar but resulted in a poor fit to the data (p- 
value=4xl0~^). 

The derivation of the A^h from the pulsar is sensitive to the 
normalization of the thermal component, dominating at low en- 
ergies. To estimate the accuracy to which we can fix the normal- 
ization of the thermal model, we calculated the ratio of bright- 
nesses in the pulsar region over the PWN region in the 0.4-0.6 
keV energy band. The resulting ratio is 1.13+0.08 confirming 
that the assumption of homogeneous thermal emission at the 
scale of the PWN is reasonable. When the normalization of the 
thermal component in the pulsar spectral model is increased by a 
factor 1.13, the best-fit value of A^h increases to 0.71+^ |^ x 10^^ 
cm"^, within statistical errors of the value quoted in Table[l] This 
test provides an estimate of the systematic error associated with 
the normalization of the thermal emission at low energies. 

To test for a possible spectral cooling in the nebula, a spec- 
trum was derived in the inner part (15"< r <60") and the outer 
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Table 1. Best-fit parameters obtained on the nebula, the pulsar and diff'erent parts o f the rim of the RX J0852. 0-4622 SNR. For each 
spectrum, a two components (absorbed VAPEC + absorbed power-law, see Sect. 2.3 for details) astrophysical background model 
was estimated from the background regions shown in Fig.[T]and|4] Only the parameters from the non-thermal components are listed 
in this Table while the parameters of the thermal component are rather constant from one region to another (kT = 0.14-0.17 keV 
and the Nw was fixed to 5 x 10^^ cm"^ for all the regions). The unabsorbed non-thermal fluxes of the PWN, PSR and SNR's rim are 
given in the 2-10 keV energy range and the A^h in units of 10^^ cm"^. The area of the extraction region and the annulus radii for the 
PSR/PWN analysis are also given. The error bars are given at 90% confidence level. 





Power-law best-fit parameters 
A^H Index Flux (10"^^ ergs/cm^/s) 


Region area 
arcmin^ 


Annulus radii 
arc sec 


Puhar &r PWN 












PSRJ0855-4644 


64+^-^^ 

^•'^^-0.11 


1 94+0.09 

^•^^-0.10 


0.26!g-03 


0.2 


0-15 


PWN 


^' ' ^-0.05 


^•^•^-0.04 




19.4 


15-150 


Inner PWN 


0.64 


-0.06 


39+0 02 

^••^^-0.06 


2.9 


15-60 


Outer PWN 


0.64 


2Q1+0.04 
^•^^-0.04 


55+o-o^ 

^-'^''^-0.06 


16.5 


60-150 


SNR 












South-East rim 




9 0+0.09 
-0.08 


1 OO+0.08 
^•"'^-0.09 


25.7 




South rim 


0.85!0-05 




o OQ+0.09 


42.1 




North-West rim 


0.68!0-04 


9 CO +0.05 
^•^'^-0.04 


4 10+0-^^ 


15.2 





part of the nebula (60"< r < 150"). Eor both regions, the as- 
trophysical background was fixed to the best-fit model obtained 
from the whole PWN in Sect. |2.3[ In order to limit the degener- 
acy between the absorption column and the photon index, the A^h 
of non-thermal emission from the nebula was fixed to the value 
derived from the pulsar. 

The best-fit parameters of the inner and outer region of the 
nebula are presented in Table [T] A significant steepening of the 
spectral index towards the external region of the nebula is ob- 
served (from 1.70+^-^^ to 2.01+^-^4) jhis is consistent with the 
scenario of a population of electrons cooling as they are advected 
away from the pulsar. 

In the view of the spectral parameter evolution in the in- 
ner nebula, the background model for the pulsar was revisited. 
The inner region of the nebula, instead of the whole nebula, was 
tested as the template background model for the pulsar. The best 
fit parameters of the pulsar with this new background model are: 
A^H=0.59!^ |4 X 10^2 cm-2 and Index=1.20!^-^^. When rescal- 
ing the normalization of the thermal component by the same 
1.13 factor (as tested in the previous section), the A^h slightly 
increases to 0.67+^ |^ x 10^^ cm'^ and Index=1.23+^ |^. Those 
parameters are compatible with the ones presented in Table [T] 
and in the following discussion only the best fit parameters from 
the pulsar presented in the table will be used. 

2.5. Spectral analysis of the rim of the SNR 

To discuss the relative distances of the SNR and of the PWN, 
we carried out a spectral analysis of diff'erent regions of the RX 
J0852.0-4622 SNR. In the same pointing where the PWN is ob- 
served, we extracted a spectrum from the rim in the region la- 
beled SE rim in Eig.[T] The region does not fully encompass the 
rim in order to avoid the bright structure of thermal emission 
from Vela that would hamper a proper determination of the ab- 
sorption column. As in the previous sections, an astrophysical 
background model (absorbed VAPEC -h absorbed power-law) is 
fitted to a region outside the SNR represented by the dashed box 
in Fig.[T] The total spectral model for the SNR's rim is the sum 



of the previously estimated background model plus an absorbed 
power-law representing the synchrotron emission from the SNR. 
The parameters of the background model are fixed in the analy- 
sis of the rim except the normalization of the thermal emission. 
The best-fit parameters of the non-thermal emission from the 
rim are listed in Table [T] The same procedure is then applied to 
other regions of the SNR where X-ray observations were avail- 
able namely the North-West and the South of the SNR. The 1.2-6 
keV image of those regions together with the spectral extraction 
regions and background regions are shown in Fig.|4] The eastern 
rim of the SNR has also been observed by XMM-Newton but the 
observation has been, unfortunately, strongly aff'ected by solar 
flares and only 5 ks out of 36 ks of observation are usable. As a 
consequence, no reliable spectral information could be derived 
from the difl'use and faint emission from the SNR in this region. 

The resulting best-fit parameters of the spectral analysis from 
three diff'erent regions of the SNR are listed in Table [T] 

2.6. Radial profiles of the SNR's rim 

To estimate the non-thermal contribution of RX J0852. 0-4622 
at the PWN position, a radial profile was extracted from a region 
of the same width and direction as the region labeled SE rim 
in Fig. [T] (right panel) but with a diff'erent length. The resulting 
profile, shown in Eig.[5](top), demonstrates that the emission of 
RX J0852.0-4622 is not negfigible at the reported distance of 
the PWN behind the shock front and its implementation in the 
astrophysical background modeling, in Sect. |2.3[ is justified. 

In order to measure the scale width of this filament, the pro- 
file was fitted with an exponentially decreasing model of the 
form ^"l^o-'^l/w^^ where d denotes the downstream side of the shock 
(same geometrical model as in Bamba et aLj |2005[ where the 



North-West rim of the SNR is studied). The best-fit model is 
overlaid in Fig. [5] (top) and the best-fit parameters is w^= 6.85 
± 0.55 arcmin. For comparison, the radial profile of the North- 
western region is shown in Fig. [5] (bottom). The latter profile was 
obtained from a region of the same direction than the white solid 
region shown in Fig. |4] (top). The same geometrical model was 
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Fig. 2. Top: Zoom on the PWN region in the 1.2-6 keV image shown 
in Fig. [T] (right). North is up and East to the left. The region is annotated 
with the annulus regions of radii Rl, R2 and R3 (15", 60" and 150" re- 
spectively) used in the morphological and spectral analysis. A contour 
at the level of 2.5 x 10~^ ph/cm^/s/arcmin^ is overlaid in red to outline 
the morphology of the inner nebula. Middle: radial profile of the PWN 
in the high energy band (1.2-6 keV) extracted from a region centered 
on the pulsar. The best-fit two components model (point-like + diffiise) 
is overlaid. Bottom: Azimuthal profile extracted from the annulus re- 
gion (15"< r <60") illustrated in Fig. [T] (right) in the high energy band. 
Position angle 0° corresponds to the West and 90° to the North. Two 
significant enhancements over the mean value (dashed line) are seen at 
~50° and ~230°. The fact that these knots of emission are separated by 
180° is suggestive of a jet morphology. 



used and the best-fit parameters are 1.61 ± 0.05 arcmin. We 
note that the downstream scale of the rim in this region is nearly 
four times thinner than in the South-East region. The implication 
of those diflferent widths on the value of the magnetic field that 
can be derived is discussed in more details in Sect. 14.51 
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Fig. 3. Top: X-ray spectrum of the nebula extracted from a 15-150" 
annulus region shown in Fig. [T| The astrophysical background com- 
ponents (dominated by: thermal emission from Vela + non- thermal 
emission from RX J0852. 0-4622) are depicted by dashed lines while 
the non-thermal emission of the nebula is illustrated by the solid line. 
Bottom: X-ray spectrum of PSR J0855-4644 extracted from a 15" ra- 
dius region. The astrophysical background (dashed lines) is composed 
of Vela's thermal emission at low energies while the contribution of the 
nebula and of RX J0852.0-4622 at higher energies is not negligible. For 
graphical purposes only the MOSl and M0S2 data are shown although 
the PN data were also used in the fitting process. 



3. Distance estimate 

The method to tackle the issue of the distance to the pulsar is 
based on the existence of a strong contrast of integrated ^^CO 
between the eastern and the western parts of the SNR, as shown 
in Fig. [7] (left). This contra st is due to the Ve la molecular ridge 



(VMR) cloud C (see Murphy & May 1991 for cloud nomen- 



clature), lying at a distance of 700+200 pc. One of the method 
to derive the latter distance is an infrared photometric study of 
the star- forming regions in the VMR (see |Liseau et al. 1992 



and references therein). Other independent methods, discussec 
in M urphy & May| ( |1991| ), give a similar distance. This cloud 
is the main contributor of ^^CO emission along the line of sight 
as seen in Fig. [6] (bottom panel) and is located at a radial veloc- 
ity in the local standard of rest of Vlsr ~ 5 km/s. If the SNR 
RX J0852.0-4622 lies behind the VMR, the observed gradient 
in ^^CO should also be visible in the measurement of the X-ray 
absorption column. 

To investigate such a correlation, absorption column density 
maps (presented in Fig.|7]) are calculated using HI data from the 



Parkes Galactic all sky survey (McClure-GriflSth s et al. 
Kalberla etZ|[2Q10| ) and ^^CO data from thcLDame et al. 



20091 
(|200T| 
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Fig. 4. XMM-Newton background subtracted flux images (combining 
MOS and PN cameras) of the North-West (Top) and South (Bottom) 
pointings of RX J0852. 0-4622 in the 1.2-6 keV energy band. As in 
Fig.[T] the solid lines illustrate the regions of spectral extraction whereas 
the dashed lines represent the regions used to estimate the astrophysi- 
cal background. The images have been adaptively smoothed in order to 
reach a S/N of 10, the color scale is in square root and the units are 
ph/cm^/s/arcmin^ . 



survey. To fully enclose the gas material lying between the ob- 
server and the back of the VMR, the data were integrated on 
radial velocities from -10 km/s to 15 km/s. Due to the intrin- 
sic velocity dispersion of the ^^CO emission line, it is unfortu- 
nately not possible to estimate the column density between the 
observer and the foreground of the VMR. The CO-to-H2 mass 
conversion factor and the HI brightness temperature to A^h used 
to produce the maps a re respectively of 1 .8x 10^^ cm"^ km~^ 
s dDame et al.l|200l|) and 1.82 x 10^^ cm'^ K"^ km"^ s jPickey 



|& Lockman[ 1 1 99 0^ The same conversion factors were used to 
compute the cumulative absorption column at the pulsar's posi- 
tion shown in Fig. [6] (top panel). From this curve, we estimate 
that, along the line of sight of the pulsar, the A^h between the 
observer and the back of the VMR is 0.97 xlO^^ cm"^ while the 
total A^H along the line of sight is 1.65 xlO^^ cm"^. 

Using the previously generated maps, a comparison plot is 
built (see Fig. [5]) where the A^h derived from the gas maps and 
integrated up to the VMR is plotted against the best-fit values 
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Fig. 5. Radial profiles obtained in the rims of the RX J0852.0-4622 
SNR in the South-East (top) and the North-West (bottom) along the 
spectral extraction regions shown in Fig. |1|4| . In the South-East, the 
distance of the PWN from the front shock is represented by the dashed 
lines. The best-fit exponentially decreasing model in the downstream 
side of the shock is illustrated with the solid lines. 



of the X-ray absorption for several positions in the sky. No cor- 
relation is observed in such a plot indicating that the SNR is 
lying in the foreground of the VMR. In addition we note that 
the A^H towards the pulsar (0.64+^ |^ x 10^^ cm"^) is slightly 
less than the column towards the SNR's rim in the same region 
(0.88^Q Q^ X 10^^ cm"^) suggesting that the pulsar is closer than 
or at the same distance as the SNR RX J0852.0-4622. 



An additional test to determine whether the pulsar is in the 
foreground or background of the VMR is to look for the ra- 
tio of neutral hydrogen (related to the X-ray absorption col- 
umn A^h) to free electrons (given by the pulsar radio dispersion 
measure DM). As discussed in Gae nsler et al.| ( [20041 ), this ratio 
has been observed to be high for pulsars which are known to 
lie behind dense molecular material (e.g. Nu/DM=^5 (40) for 
PSRJ1747-2958 and PSR 1757-24 respectively) whereas the 
typical value for other pulsars is of the order of 5-10. The value 
we derive for PSR J0855-4644 is 8+1.5 which suggests that the 
pulsar is not located behind dense material. 

We therefore conclude that both the SNR and the pulsar are 
lying in the foreground of the VMR i.e. at a distance d < 900 
pc. This distance constraint to the RX J0852.0-4622 SNR is 
compatible with a recent estimate based on the measurement of 
the proper motion of the shock in X-rays in the North-western 
part of the rim by |Katsuda et al.| ( [2008| ) who derived d~750 pc. 
Concerning the pulsar, a closer distance (than the 4 kpc derived 
from the radio DM measurements) confirms a suggestion made 
by [Redman & Meaburn] ( [2005] ) who derived a distance of -750 
pc. 
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□h [10'^ cm-'] 

Fig. 7. Maps of the absorption column density, derived from the molecular gas (left panel), the atomic gas (middle panel) and the sum of the two 
previous components (right panel) when integrating over radial velocities from -10 km/s to 15 km/s (see Sect.|3]for details). This specific velocity 
range provides an estimate of the absorption column for a source located just behind the Vela molecular ridge. The position of PSR J0855-4644 
is shown with a black cross and the XMM-Newton pointings used in this study are represented by the white circles in the right panel. The white 
contours are taken from a ROSAT high energy image (E>1.3 keV) and reveal the SNR RX J0852.0-4622. 




derived from the X-ray analysis show high values ranging from 
0.64 to 0.88x10^^ cm"^. This is intriguing as the Nu towards the 



50 100 
Velocity (LSR) km/s 
Fig. 6. Top: cumulative absorption column density (dot-dashed 
line) together with the ^^CO and HI contributions as a function 
of the radial velocity integrated on a 0.1° box region centered 
on the pulsar. Bottom: corresponding ^^CO and HI spectra at 
the position of the pulsar. The Vela molecular ridge, centered at 
Vlsr~5 km/s, is the main ^^CO structure along the pulsar's line 
of sight. 



4. Discussion 

4.1. Proposed scenario 

In the previous sections, we have shown that RX J0852. 0-4622 
and PSR J0855-4644 are lying in the foreground or at the inter- 
face of the VMR. The SNR and the pulsar are therefore thought 
to be nearby (d< 900 pc) but surprisingly the absorption columns 



Vela SNR whose distance is known precisel}rlto be d=287^ 



19 , 



( Dodson et al. ,12003') is much lower at ~10 



cm 



To understand this apparent discrepancy, we investigated the 
spatial distribution and the density of the local gas towards the 
Vela region. We note that due to the intrinsic velocity dispersion 
of the VMR emission both in HI and in ^^CO (see Fig. [6] bottom) 
and the proximity of the object, it is not possible to disentan- 
gle the emission from the VMR from the foreground emission. 
However, the density of the nearby gas (the Local Bubble) can 
be probed by other methods such as the study of the absorption 
line fe atures of Nal and Call towards stars with known distances 
(see Lalleme nt et"aLl [20031 [Welsh et al.l|201 0). These studies 
revealed interstellar tunnels of low density connecting the Local 
Bubble to other surroundings cavities. One of these tunnels is 
in the line of s ight of the Vela r egion at ^ ~ 265° extending at 
least to 250 pc ([Welsh et al.||2010 ). This tunnel could explain the 
low value of the A^h measured in X-rays in the direction of the 
Vela SNR. Using the A^h measured towards the Vela SNR (10^^ 
cm"^), the mean density in this line of sight is of the order of 0.1 
cm"^. 

After the Local Bubble, the line of sight encounters the Gum 
nebula, a large scale structure centered at f = 262°, b = -3° 
and located at a distance of ~ 500 pc for a physical radius of 
~ 150 pc and an angular radius of 17°. Over a distance of 200 or 
600 pc (distance between the Vela SNR at 287 pc and the VMR 
at 700+200 pc), the required densities must be of the order of 
9 to 3 cm"^ respectively to reach the A^h value measured in X- 
rays of ~ 0.6 x 10^^ cm'^ towards PSR J0855-4644. When inte- 
grating the HI spectrum between the Earth and the background 
of the VMR ( -10 km/s to 15 km/s), the HI alone contributes 
to 0.3 X 10^^ cm"^ in the line of sight towards the pulsar cor- 
responding to a mean density of 2-5 cm"^ over a distance of 
600 and 200 pc respectively. In addition the HI emission is ho- 
mogeneously spatially distributed around the RX J0852. 0-4622 



The distance to the Vela pulsar has been measured by parallax. 
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SNR region and is not correlated with the VMR as shown in the 
HI map shown in Fig. [7] (middle panel) thus contributing to all 
the lines of sight used in the X-ray analysis which is consistent 
with the relatively high A^h values measured across the SNR. The 
densities aforementioned are compatible with the density of ~7 
cm~^ measured in the inner region of the Gum nebula (the thick 



HI shell discussed in Dubner et al. 1992[ Reynoso & Dubner 



|1997) that is spatially overlapping our region of interest. Hq, ob- 



servations of the nebula by [Reynolds] ( | 1 976| ) provide similar val- 
ues for the densities (1 to 5 cm""^). One possible origin of the 
Gum nebula is that the supershell is produced by the combina- 
tion of stellar winds and r epeated SN explosions (McC ray ~&\ 
Kafatos 1987[ Reynoso & Dubner. 1991) and it is therefore a 



plausible location for PSR J0855-4644 and RX J0852.0-4622. 

Finally, changing the electron density of the Gum nebula, in 
the NE2001 model, from 0.5 cm"^ (the currently implemented 
uniform value for the Gum nebula) to 1 cm"^ is sufficient to 
reconcile a DM of 239 cm"^ pc with a distance to the pulsar 
of 900 pc. This does not seem unreasonable with regards to the 
possible large density variations in the Gum nebula. 

In summary, the line of sight towards the pulsar encounters 
first a tunnel of low density gas (-0.1 cm"^) up to the Vela SNR 
then enters the Gum nebula thick HI shell with densities of the 
order of 1-5 cm"^. This succession of low and high density en- 
vironment can explain the low A^h towards the Vela SNR and 
the high A^h towards the RX J0852.0-4622 SNR and the pulsar. 
While the HI alone contributes to 0.3 x 10^^ cm"^, additional 
molecular material in the line of sight could contribute to the 
measured X-ray absorption either in the foreground of the VMR[^ 
or from the external regions of the VMR if the pulsar and SNR 
are lying at the interface of the molecular cloud. 



4.2. Is the pulsar associated with the RX J0852.0-4622 
SNR ? 

Although the pulsar and the SNR may be at similar distances, an 
association is unlikely due to the characteristic age of the pul- 
sar (Tc=140 kyrs) being much larger than the SNR age (tsNR ~ 
4000 yrs). Assuming that could be very diff'erent from the 
real age, the required pulsar's kick velocity to explain its cur- 
rent position would be -3000 km/s in order to travel 12 pc from 
the SNR center (1° at 750 pc) for a SNR age of 400 years. Th is 
value is much higher than observed in other pulsars ("Hobbs et al. , 
[2005 ) . Moreover for such a high shock speed a bow-shock or a 
tail structure should appear in the PWN which is not the case. 
In addition, a candidate neutron star CXOU J085201. 4-461753 
2002| ) lies near the geometrical center 



(see e.g.: Kargaltsev et al. 



of the SNR RX J0852.0-4622 which appears as a more simple 
scenario. 

4.3. TeV emission from the PWN 

With the new distance derived in this study, the ratio E/d^ for 
PSR J0855-4644 reaches a value > 10^^ erg/s/kpc^ (for d< 900 
pc and E = 1.1 x 10^^ ergs). According to the TeV y-ray PWN 
population study presented by Carrigan et al. ( 2008 ), such a high 
E/d'^ pulsar should be detectable by Cherenkov telescopes in 
y-rays. The region of the RX J0852.0-4622 SNR has been ex- 
tensively studied by the H.E.S.S. experiment and an excess of 
TeV y-rays is observed at the position of the pulsar (see Fig. 3 



Paz Arribas et al. 2012). However as the pulsar is in spatial 



coincidence with the shell of the SNR, it is not possible, given 
the spatial resolution of the H.E.S.S. telescopes, to disentangle 
the contribution from the shell of the SNR from the contribu- 
tion of the PWN. In addition, the azimuthal profile of the y-ray 
emission along the shell of the SNR does not reveal any sig- 
nificant emission increase with regards to the surrounding SNR 
shell emission at the position of pulsar ( Paz Arribas et al. |2012 ). 



4.4. Pulsar contribution to observed e /e+ CRs 

This new distance estimate of the pulsar strongly diff'ers from 
the distance estimated with the radio dispersion measure (d=4 
kpc) and therefore implies that PSR J0855 -4644 is an energetic 
and nearby pulsar. After the Vela pulsar, it is the most energetic 
pulsar within 1 kpc and could therefore contribute to the spec- 



trum of cosmic-rays e /e^ received at Earth (see Delahaye et al. 



2010, for the contribution of nearby pulsars). Obtaining a gooc 
knowledge of the nearby pulsar population is a necessary step to 
understand the rise in positron fraction with increasin g energy 
measured by the PAMELA experiment (Adriani et al., 2009) 
and confirmed by the Fermi space telescope ( [Ackermann et al.] 
,2012 ). 

As the pulsar is rather old (Tc= 140 kyrs) and no traces 
of the host SNR have been observed, one can reasonably as- 
sume that the CR accelerated in the PWN have had time to 
escape the SNR. The diff'usion timescale of those escaping 
CRs to reach the Earth scales as tp^opiE) ~ d^/D(E) where 
D(E) is the diff'usion coefficient with a typical value of ~ 
3 X 10^^(£^/1 GeV)^ cm^s with a large uncertainty on the en- 
ergy dependence ^=0.3-0.6 (see [Strong et al.[ |2007[ for a re- 
view of CR propagation in the Galaxy). Therefore, tpYop(E) ~ 
100 kyrs ((i/100pc)2(^/lGeV)-^ ^ For a typical value of ^=0.5, 
neither PSR J0855-4644 nor the Vela pulsar (Tc= 11 kyrs) are 
thought to contribute to the observed flux of e"/e^ at 10 GeV 
as the difl'usion timescale is respectively 1500 and 280 kyrs (for 
d=750 and 290 pc respectively). By equating ^prop = Tc and for 
^=0.5, the energy of the particles that could have difl'used up to 
Earth is respectively 1 TeV and 6 TeV for PSR J0855-4644 and 
the Vela pulsar. 

The previous estimates assume that the particles start to dif- 
fuse at the birth of the SNR and do not take into account the 
time for the particles accelerated in the PWN to break through 
the shell of the SNR. This escape effectcould be not negligible 



in the case of the Vela pulsar (see Hinton et al. 2011, for a study 



^ Such material can unfortunately not be disentangled from the VMR 
in the ^^CO spectrum. 



of the difl'usive escape efl'ect on the particle spectrum) due to the 
relative young age of the system (tc= 11 kyrs). 



4.5. Comparison of the width of the rims of the SNR 

In the XMM-Newton observation of the South-East region of the 
RX J0852. 0-4622 SNR, a faint rim of non-thermal emission was 
detected (see Fig. [T] right). In Sect. 2^ the downstream thick- 
ness of the rim was measured with an exponentially decreasing 
model. We observed that the filament in the South-East is nearly 
four times thicker (wd= 6.85 ± 0.55 arcmin) than the bright and 
thin (wd= 1.61 ± 0.05 arcmin) filament observed in the North- 
West of the SNR shown in Fig. |4] 

Assuming that the thickness of the rim is limited by the syn- 
chrotron losses while the electrons are advected away from the 
shock, one can derive a lower limit to the value of the down- 
stream magnetic field B^. To derive such a value, we adapted 
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Fig. 8. Comparison plot of the absorption column density values de- 
rived from the atomic and molecular gas and from the X-rays using 
XMM -Newton observations. The X-ray columns derived towards the 
rims of the RX J0852.0-4622 SNR, PSR J0855-4644 and the PWN are 
represented by black circles, a blue square and a green square respec- 
tively. The errors on the X-ray A^h are 1 cr error bars. The label over 
each point corresponds to regions where the values have been measured 
as defined in Fig. [7] 



the equation 9 from jParizot et al.| ( |2QQ6| ) to keep track of the 
observed photon energy E^h and obtained the following relation: 



ph,keV 



(1) 



where ysh,3 is the shock speed in units of 1000 km/s, 
^obs,arcmin is the obscrvcd width of the rim at half maximum 
in units of arc minutes (^obs=>^d x log(2)) and £^ph the aver- 
age emitted photon energy in keV. For a distance to the RX 
J0852.0-4622 SNR of 750 pc (see 



Katsuda et al. 2008 and 



Sect. |3] of this article for a more detailed dis cussion of the S NR 
distance), a shock speed of 3000 km/s ( Katsuda et al.||2008|) and 
a mean photon energy of ~2.5 keV, the downstream magnetic 
field value is 24 yuG and 63 yuG respectively for the South-East 
and for the North-West rim. 

We note that in the North- West region, an X-ray proper mo- 
tion of the shock (Osh=Vsh/d) has been measureqjby Katsuda 
et al.| ( |2QQ8| . As the ratio of V^h/d can be replaced by 4h in Eq. 
(1), the estimate of the magnetic field, at least in the North-West, 
is in fact independent of the hypothesis on the shock speed and 
distance. 

The magnetic field downstream of the shock is therefore sig- 
nificantly diflTerent in those two regions. This emphasizes the fact 
that the magnetic field derived from a one-zone modeling of the 
SED of the whole SNR (giving an average value of the magnetic 
field) is not directly comparable to the magnetic field derived 
from a filament in a specific region of the SNR. 

5. Conclusions 

The search for the X-ray counterpart of PSR J0855-4644 using 
XMM-Newton led to the following conclusions: 

1 . A point-like X-ray counterpart has been discovered in spa- 
tial coincidence with the position of the radio pulsar and a 
surrounding extended emission of non-thermal origin was 



^ The value of Vsh and d used in the previous calculation were derived 
from the aforementioned proper motion measurements. 



found. The nebula has a maximum extent of the order of 
150" and the azimuthal profile has revealed two knots of 
emission (with and angular separation of -180°), suggest- 
ing possible jets around the pulsar. High spatial resolution 
observations with Chandra might reveal the fine structure of 
the nebula. 

2. A faint rim of non-thermal emission representing the South- 
Eastem part of the RX J0852.0-4622 SNR was detected. 
The filament in this region is thicker than the bright and 
thin X-ray filaments observed in the North-West of the SNR. 
Assuming that the thickness of the rim is dominated by the 
synchrotron losses, the downstream magnetic field is 24 yuG 
and 63 yuG respectively for the South-East and for the North- 
West rim. Those significantly diflTerent values might indicate 
diflTerent acceleration and/or environmental conditions at the 
shock around the SNR. 

3. By comparing the X-ray column density on the pulsar and on 
diflTerent regions of the SNR with the gas map (using ^^CO 
and HI emission as tracers), we concluded that both the RX 
J0852.0-4622 SNR and PSR J0855-4644 are in the fore- 
ground or at the interface of the Vela Molecular Ridge re- 
sulting in an upper limit of d< 900 pc. Although the pulsar 
and the SNR may be at similar distances, they are not likely 
to be associated. 

4. The new estimate of the pulsar distance completes our 
knowledge of the nearby pulsar population. If the distance 
of the pulsar is confirmed through e.g. radio parallax mea- 
surements, PSR J0855-4644 would be the second most en- 
ergetic pulsar (after the Vela pulsar) within 1 kpc and may 
contribute to the observed e"/e^ spectrum. 
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Appendix A: Calibration issues with the PN camera 

The thermal emission from the Vela SNR is extremely bright at 
low energies (E<1 keV), an energy range where the PN camera 
is more sensitive than the MOS. With such a statistics we are 
limited by our knowledge of the calibration of the camera. Here 
we have observed a deficiency in the low energy calibration of 
the spectral response of the PN camera. This eflfect is most visi- 
ble around 0.6 keV in Fig. |A.1| where the Oxygen lines can not 
be reproduced by the APEC thermal model. 

This is a known calibration issu ^that is being investigated 
by the XMM-Newton PN calibration team. There are no known 
similar calibration issues for the MOS instruments. Following 
their recommendation, we allowed for an empirical shift in en- 
ergy in the PN response file using the gain response model im- 
plemented in Xspec. The fitted energy shift is of the order of 15 
eV. Similar best-fit parameters where obtained when fitting in- 
dependently the MOS and response corrected PN spectra. This 
confirms that the correction of the response matrix does not alter 
the physical parameters that we derive from the spectrum. In the 
spectral analysis of the article we have therefore used the MOS 1, 
M0S2 and response corrected PN spectra. 



^ More information on this calibration issue can be found here: 
http : //w w w2 .le . ac . uk/departments/phy sic s/res earch/ src/Mis sions/xmm- 1 
,newton/technical/cal-meetings jimadrid-03 1 0, 
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Energy (keV) ^ 

Fig. A.l. Top: X-ray spectrum of the nebula extracted from a 15-150" 
annulus region shown in Fig. [2] focused on the thermal emission with 
MOSl, M0S2 and PN spectra represented in black, red and green re- 
spectively. The same spectral model as in Fig. [3] is overlaid. Bottom: 
Same spectra as in the Top panel but with the PN response matrix mod- 
ified by the gain response model implemented in Xspec. 
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